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Electronic and magnetic structures of the newly synthesized cobalt oxyselenide La2Co2Se2 03 
(structurally similar to the superconducting iron pnictides) are studied through density functional 
calculations. The obtained results show that this material is a Mott insulator, and that it has a very 
stable Co^^ high-spin ground state with a i29-like orbital ordering, which is substantiated by the 
calculated crystal-field excitation energies. The square lattice of the Co"^^ spins is found to have a 
strong antiferro (a weak ferro) magnetic coupling for the second nearest neighbors (2nn) via O (862) 
and an intermediate antiferro one for the Inn, with the strength ratio about 10:1:3. The present 
results account for the available experimental data of magnetism, and the prediction of a planar 
frustrated (2x2) antiferromagnetic structure would motivate a new experiment. 

PACS numbers: 75.30.-m, 71.20.-b, 71.27.-|-a, 71.70.-d 
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Very recently transition-metal oxychalcogenides re- 
ceive revived interest^"— due to their structural similar- 
ity to the superconducting iron pnictides^ in terms of 
a common square lattice of the transition-metal species 
[with a local octahedral (tetrahedral) coordination in the 
former (latter)]. In contrast to bad metal properties of 
the parent iron pnictides, a Mott insulating behavior has 
been found in the iron or cobalt oxychalcogenides i^i^ 
As the strength of electron correlations in the pnictides 
is currently under a hot debate, and as the supercon- 
ductivity and magnetism therein have a close relation- 
ship, it is worth a lot to study promptly the correlated- 
electron structure and magnetic properties of those Mott 
insulating oxychalcogenides, which would help to under- 
stand the pnictides upon their approaching the limit of a 
Mott localization. The structure and magnetic properties 
of (Ba,Sr)2Fe2(Se,S)20F2 have recently been studied by 
Kabbour et al.^ The observed long-range antiferromag- 
netic (AF) ordering below 83.6-106.2 K was proposed to 
be a rare example of a frustrated AF checkerboard spin 
lattice. More recently, Zhu et al. studied a Mott localiza- 
tion in La2Fe2(Se, 8)203 and found that those materials 
have moderate charge gaps.— Their results support the 
notion that the iron pnictides and chalcogenides possess 
intermediately strong electron correlations and are not 
too far away from Mott localizations. 

In this Rapid Communication, we study the cobalt 
oxyselenide La2Co2Se203, which was newly synthesized 
by Wang et al and found by them to be insulating and 
to have an AF transition at Tap ^220 K.^ Those au- 
thors suggested an unusual low-spin (LS, 5=1/2) state 
of the Co^"*" ions and a corresponding orbital order for 
understanding of the insulating behavior and AF in the 
checkerboard spin-lattice. The Co^+ LS state is quite 
surprising, in view of the common high-spin (HS, S'=3/2) 
state. Having a close look at the calculated results of 
WangeiaZ. (Fig. 6 in Ref . 4^ using the local-spin-density 
approximation plus Hubbard U (LSDA-f-C/) method, one 
may wonder why three singly occupied orthogonal or- 
bitals have two up-spins and one down-spin, which ob- 



viously violates the Hund's first rule. This infers that 
their LSDA+U calculation with an orbital-polarized po- 
tential ran into an exotic solution which should not be 
the ground state. Secondly, if the Co^"*" is in the LS 
state as Wang et al concluded using measurements of 
magnetic susceptibility and specific heat^ (which will be 
commented on below), the magnetic coupling between 
the S'=l/2 ions, with a most probable frustration, would 
not be expected to yield the quite high Tap ~220 K. 
Thus the electronic structure and magnetism of this in- 
teresting new material were poorly understood. For this 
reason, we carry out a rather complete set of density func- 
tional calculations to study the electronic and magnetic 
structures of La2Co2Se203, which are closely related to 
the Co^"*" spin and orbital states. In strong contrast to 
the previous finding^ all the present calculations using 
LSDA, LSDA-Fock-0.25, and LSBA+U functionals con- 
sistently conclude the HS ground state, which is also sub- 
stantiated by the calculated crystal-field level splittings. 
The orbital multiplet effect and a t2g-hke orbital ordering 
arc discussed. Moreover, the exchange constants of the 
Co^+ square lattice are calculated and a planar frustrated 
(2x2) AF structure is proposed. The present work sheds 
new light on this interesting material, and a comparison 
is also made with the above iron oxychalcogenides. 

La2Co2Se203 has the Co2Se20 layers which are sepa- 
rated by the La202 layers, see Fig. 1(a). Each Co2Se20 
layer consists of the CoSe402 octahedra having the Co-0 
(Co-Se) bondlength of 2.035 (2.688) AA Those octahedra 
are face-shared for the first-nearest-neighboring Co sites 
and edge- (corner-) shared via Se2 (O) for the second- 
nearest Co neighbors. In our calculations, we used for 
each CoSe402 octahedron such corresponding local xyz 
coordinates that for the body-centered Co site as an ex- 
ample, the z axis is along the Co-0 bond and the x along 
the crystallographic c axis. As a result of the D2h point- 
group symmetry seen by the Co^"*" ions, the lower- lying 
t2g-like orbitals {x'^-y^, yz, and xz) and the higher eg [xy 
and iz^-r"^) all split,— see Fig. 1(b) and the detailed re- 
sults below. Our calculations were carried out using the 
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FIG. 1: (Color online) (a) Crystal structure of La2Co2Se203 
having the Co2Se20 layers. For each constituent CoSe402 oc- 
tahedron the corresponding local coordinates are used, e.g., 
the xyz for the body-centered Co site and the y and 2 inter- 
changed for the c-axis edge Co^^. (b) Orbitally resolved Co 
3d DOS for the nonmagnetic state by LDA. The lower-lying 
i29-like orbitals (x^-y^, yz, and xz) and the higher eg {xy and 
3z^-r^) all split, (c) DOS of the HS FM half-metallic state 
by LSDA. The solid (dashed) curves refer to the majority 
(minority) spin. Fermi level is set at zero energy. 



full-potential augmented plane waves plus local orbital 
method.* The muffin-tin sphere radii are chosen to be 
2.8, 2.2, 2.0 and 1.5 Bohr for La, Co, Se, and O atoms, 
respectively. The cut-off energy of 16 Ryd is used for 
plane wave expansion, and 15x15x3 k mesh for integra- 
tions over the Brillouin zone, both of which ensure a high 
accuracy of our results. 

We start with the standard LSDA calculations. The 
results show that independent of the initialized spin mo- 
ment, a ferromagnetic (FM) state always converges to 
the same HS solution, see its density of states (DOS) 
in Fig. 1(c). It is half metallic, with the fully occupied 
(and thus insulating) spin-majority channel and partially 
occupied (metallic) Co M bands. This solution is calcu- 
lated to be more stable than the nonmagnetic state by 
593 meV/f.u. We also did for the assumed LS state a 
fixed-spin-moment calculation (otherwise it converges to 
the above HS ground state) and find that the LS-fixed 
state is less stable than the HS state by 435 meV/f.u., see 
Table I. Moreover, our calculations show that the HS FM 
state is more stable than the HS checkerboard- AF state 
by 26 meV/f.u. This is most probably due to a max- 
imal kinetic energy gain in the FM half-metallic state. 
Note, however, that this FM half metallic ground state 
disagrees with the experimental insulating and AF state. 



TABLE I: Relative total energies in unit of meV/f.u. (and the 
Co spin moment in unit of /is) of La2Co2 86203 having the 
HS or LS Co^+ ions in the checkerboard AF or FM state. The 
three functionals all show the HS ground state. The symbol 
"— >" stands for "converges to". 
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258 (2.65) 
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By a comparison with LaFeAsO (La2Fe2Se203) in which 
the Fe 3d bandwidth is 4.4 (3.2) eV,^. the entire Co 3d 
bandwidth in La2Co2Se203 is only 2.8 eV, see Fig. 1(b). 
This signals a stronger correlation effect in La2Co2Se203 
due to the band narrowing.— Indeed a presence of the 
local magnetic moment suggests a Mott localization be- 
havior of La2Co2Se203.'* Therefore, it is not surprising 
that the standard LSDA calculations (albeit quite infor- 
mative) cannot provide a correct description of the elec- 
tronic and magnetic structure of this Mott insulator. 

In order to study the electron correlation effect, we car- 
ried out hybrid functional calculations using the LSDA- 
Fock-0.25 which mixes 25% Fock exact exchange into the 
LSDA exchange functional and is applied to the Co 3d 
electrons. This was demonstrated to be an alternative 
way (to the LSDA±[/, see below) for studying correlated 
electron systems.— »i^ Our hybrid functional calculations 
show again that the Co^"*" HS is the ground state, be- 
ing much more stable than the LS state by about 1.5 
eV/f.u., see Table I. [The corresponding much smaller 
value of about 0.4 eV/f.u. given by the LSDA may well 
be due to the artificial (half-) metallic band structure 
which significantly suppresses the energy splitting of the 
spin and orbital multiplets of the Co^"^ ions.—] Fig. 2 
shows the orbitally resolved DOS for the HS (left panel) 
and LS (right panel) solutions in the checkerboard AF 
state. The HS Mott insulating solution has one t2g hole 
on the minority-spin xz orbital and two Cg {xy and 3z^- 
r^) holes, thus forming a t2g orbital ordering. This is the 
ground-state solution of the orbital multiplets for the HS 
Co^"*" ions, as confirmed by the following LSDA+[/ calcu- 
lations. In this HS insulating state, the AF exchange and 
superexchange are expected to be dominant, and indeed 
our calculations show that the HS checkerboard AF state 
is more stable than the FM state by 258 meV/f.u., see 
Table I. In contrast, the very unstable LS state has the 
^2q^q configuration with the e]. being an xy electron, fully 
according with the crystal-field level sequence calculated 
below. The LS AF state is (a little) more stable than the 
LS FM state, but the energy difference of 24 meV/f.u. is 
less than one tenth of that for the HS state. Although ex- 
change constants may well be overestimated by a hybrid 
functionalj= those two values confirm, as expected, that 
the LS state would yield much weaker magnetic couplings 
than the HS state. 
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FIG. 2: DOS of the HS (left panel) and LS (right panel) 
solutions in the checkerboard AF state by the LSDA-Fock- 
0.25 hybrid functional. The Co^^ ions in the HS ground state 
have an xz-type (f2g hole) orbital ordering. The much less 
stable LS state has the t2g+e]i{:xy) configuration. See Fig. 1 
for other notations. 



As the LSDA-fC/ method^'^ has been widely used to 
study the electronic and magnetic properties of the cor- 
related electron systems, here we further performed a set 
of LSDA+C/ calculations using the double-counting cor- 
rection in an atomic limit and the effective C/=4 eVJ^ 
We have initialized the corresponding density matrix 
(and orbital-dependent potential) and then done self- 
consistently a full electronic relaxation. This procedure 
turns out to be suitable for study of the spin and or- 
bital physics. ^^ The obtained electronic structure, the HS 
ground state, the orbital multiplet splitting, the exchange 
constants and the ground-state magnetic structure are 
discussed as follows. 

For both the FM and checkerboard AF configurations, 
the LS state turns out to be unstable and converges 
to the HS ground state. For the HS ground state, we 
tested all the ten different orbital configurations of the 
two minority-spin electrons (the majority-spin channel is 
closed in the Co^+ HS ground state). We find that the 
five lowest-energy configurations states, the {x'^-y'^){yz), 
ix'^-y'^){xz), iyz){xz), {yz){xy), and {x'^-y'^){^z'^-r'^), or- 
der as 0, 60, 457, 750, and 790 meV/Co (see Fig. 1 for 
orbital indices). Thus, the ground state of the orbital 
multiplets is t\g^^e'^g^{x'^-y'^Yi{yz)\, confirming the above 
hybrid functional results. The corresponding DOS re- 
sults (not shown here) are very similar to the latter. As 
each Co'^^ ion has a t2g hole on the xz orbital in its 
local coordinates and the local coordinates are rotated 
around the crystallographic c-axis by 7r/2, to say, from 
the body-centered Co^+ site to the c-axis edge one [see 
Fig. 1(a)], each Co2Se20 layer has an antiferro a;z-type 
hole-orbital ordering. Moreover, from those above en- 
ergy differences and by neglecting the difference of the 
inter-orbital Coulomb interactions, the crystal-field level 
splittings (or to say, crystal-field excitation energies) are 
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J\ Ji Jz 

U=4 -9.97 -3.06 0.98 

^=3 -12.39 -4.04 0.84 
L/=2 -16.32 -5.36 0.78 



FIG. 3: (Color online) The ground-state (2x2) planar mag- 
netic structure of the HS Co^^ square lattice, with the 2nn 
AF exchange Ji via O (red bold lines), the Inn AF J2 (blue 
solid lines), and the 2nn FM exchange J3 via Se2 (cyan thin 
lines). See also Fig. 1(a). The J2 couplings are frustrated. 
Only the spin-up and -down Co'^"'" sites are shown for clarity. 
The listed exchange constants (in unit of meV) are calculated 
using the LSDA-f (/ with effective f/=4, 3, or 2 eV. 



estimated: the x'^-y^, yz, xz, xy, and 3z^-r^ order as 
0, 397, 457, 750, and 1187 meV. By a comparison with 
the LS state [t^gixyY), the HS ground state ends with 
an xz -^'iz'^-r^ transition which gains, by a simple es- 
timate, two times Hund exchange (in total about 2 eV) 
with a much lower cost of the crystal-field energy of 730 
meV. This is the reason that the HS is the very stable 
ground state of this Mott insulator, having a localized 
spin moment larger than 2.5 jjlb as seen in Table Ij^^ 

Now we estimate the three exchange constants of the 
HS Co^"*" ions in the square lattice (Fig. 3), by calculat- 
ing four different magnetic structures (the same as used 
in Ref. \^ see Fig. 9 therein) and then mapping the 
total-energy differences (i.e., exchange energies) onto a 
Heisenberg model. The resultant values are -9.97 meV 
(Ji, AF, counted once per pair), -3.06 meV (J2), and 
0.98 meV (J3, FM) for the second-nearest-neighboring 
(2nn) Co pair via O, the first-nearest-neighboring (Inn) 
pair, and the 2nn pair via Se2, respectively, as shown in 
Fig. 3. As the Ji is more than three times as big as the 
J2 (well meeting the requirement | Ji | > | J2 1 ) and the J3 
is positive, the ground state AF structure [Eex = (^1- 
Jz)S'^ per Co^"*" with S=2>/2] is what we show in Fig. 3: 
It is more stable (albeit the frustrated J2 couplings) than 
the checkerboard AF state [E^x = {-Ji+'2J2^J3)S^ per 
Co] by 31 meV/Co, and than the stripe AF state [E^x 
= {Ji+J3)S^] by 4 meV/Co. We also calculated those 
J values using a smaller U (=3 or 2 eV) when the sys- 
tem approaches an itinerant limit, and find that while 
the FM J3 becomes a bit smaller, both the AF Ji and J2 
get larger with a nearly constant Ji/ J2 ratio of about 3. 
Thus, the HS planar (2x2) AF structure as shown in Fig. 
3 remains to be the ground state as the system is a Mott 
insulator with the effective U > 2 eV. This prediction 
calls for an experimental study. 

Here we compare our LSDA-|-[/ results with those cal- 
culated previously for the isostructural Ba2Fe2Se20F2 
and La2Fe2Se203^i^ All the three materials are calcu- 
lated to be a Mott insulator with a band gap opening 



by Hubbard C/, and in the Mott insulating phase each of 
the three exchange constants has the same sign (AF or 
FM) in the three materials. Both the Fe2+ (S'Fe=2) and 
Co^"*" (5*00=3/2) ions are found to be in the HS state. 
Considering the strongest AF coupling for the Co or Fe 
pair via O, JiSq^ in La2Co2Se203 is about 1.5 times 
as big as JiSj^ in BaaFeaSeaOFa (Ref. I2!) and J^S^.^ 
in La2Fe2Se203 (Ref. 5) (using those J values at the 
U~3 eV for a direct comparison). Moreover, the mag- 
netic frustration appears weaker in La2Co2 86203 due to 
the smaller ratio J2/J1 < 1/3 (the corresponding ratio 
>0.5 in Ba2Fe2Se20F2 and La2Fc2 80203,^1^) Both fac- 
tors account for the observations that the Taf ~220 K 
in La2Co2Se203 (Ref. 0) is about two times as high 
as that of around 100 K in both Ba2Fe2Se20F2 and 
La2Fe2Se203.^'^ Note that if the 00^+ was in the LS 
state, the Tap would decrease probably by one order of 
magnitude as discussed above by the hybrid functional 
calculations. This would obviously disagree with the ex- 
periments. Therefore the LS-state model is discarded by 
the present work. 

Now we briefly comment on the recent experiments 
carried out by Wang et al. for La2 00286203.'* They used 
a Curie- Weiss law to fit the magnetic susceptibility in the 
quite narrow range of317K<T< 400 K, and they also 
estimated a change of the magnetic entropy below 255 
K, both of which led them to a conclusion of the unusual 
L8 state. As this layered material has a two-dimensional 
spin square-lattice with a considerably strong magnetic 
frustration as found in the above calculations, the ob- 
served Tap ^ 220 K may well be lower than the domi- 
nating AF coupling strength which is suppressed by other 
frustrated magnetic interactions. This brings about a 



broad maximum of the susceptibility above the Tap as 
observed, and thus the above fitting using a Curie- Weiss 
law is less reliable: the strongest AF coupling is not yet 
well broken, leading to an underestimate of the effective 
moment. For the same reason, the decreasing magnetic 
entropy below 255 K represents only a partial but not 
a full contribution of the spin degree of freedom, leading 
again to an underestimate of the spin degeneracy. As the 
H8 ground state is very stable as we find above, and it 
has no room for a transition to a higher spin state with 
increasing temperature, a magnetic susceptibility mea- 
surement at a higher temperature, e.g., up to about 600 
K, would be able to check it. 

To conclude, we studied the electronic structure, spin 
state, and magnetism of the newly synthesized square- 
lattice cobalt oxyselenide La2Co28e203, using the L8DA, 
L8DA-Fock-0.25, and L8DA-HC/ calculations. While the 
L8DA predicts this material to be a possible ferromag- 
netic half metal, both the L8DA-Fock-0.25 hybrid func- 
tional and L8DA-I-C/ calculations show a Mott insulating 
and antiferromagnetic ground state in agreement with 
experiments, and a t2g orbital ordering as well. The 
exchange constants of the frustrated antiferromagnetism 
are calculated. In contrast to the newly suggested un- 
usual low-spin state of the Co^+ ions, all the present cal- 
culations conclude the very stable high-spin ground state, 
which is substantiated by the calculated crystal-field level 
splittings. A comparison with the iron oxychalcogenides 
is made, and the experiments are explained. 

The author thanks Thomas Lorenz for helpful dis- 
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